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ABSTRACT 

Neutrino-driven wind from young hot neutron star, which is formed by supernova 
explosion, is the most promising candidate site for r-process nucleosynthesis. We study 
general relativistic effects on this wind in Schwarzschild geometry in order to look for 
suitable conditions for a successful r-process nucleosynthesis. It is quantitatively dis- 
cussed that the general relativistic effects play a significant role in increasing entropy 
and decreasing dynamic time scale of the neutrino-driven wind. Exploring wide param- 
eter region which determines the expansion dynamics of the wind, we find interesting 
physical conditions which lead to successful r-process nucleosynthesis. The conditions 
which we found realize in the neutrino-driven wind with very short dynamic time scale 
Tdyn ~ 6 ms and relatively low entropy S ~ 140. We carry out the a-process and r- 
process nucleosynthesis calculation on these conditions by the use of our single network 
code including over 3000 isotopes, and confirm quantitatively that the second and third 
r-process abundance peaks are produced in the neutrino-driven wind. 
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1. Introduction 

The r-process is a nucleosyn thesis process to pro - 
duce elements heavier than iron( Burbidge et al. 1957 ). 
They occupy nearly half of the massive nuclear species, 
and show typical abundance peaks around nuclear 
masses A=:80, 130 and 195, whose neutron numbers 
are slightly smaller than the magic numbers N=50, 82 



an d 126, respectively. This fact auggcata that the r - 



ies of neutrino-driven wind followed the successful de- 
tection of energetic neutrinos from SN1987A, which 
raised the possibility of finding the r-process nucle- 
osynthesis in this wind. 

Although there are several numerical simulations 
of the neutrino-driven wind, results are very different 
from one another, depending on models and meth 
ods adopted in literature (Woosley et al. 1994, Witti 



Janka, & Takahashi 1994, Takahashi, Witti, fc Janka] 



pro CC33 elements have completely different origin from 
the^ 3 - procc33 clcmcnta whose abundance peaks arc lo - 
cated just on the neutron magic numbers. The r- 
process elements are presumed to be produced in an 



exp loaivc environment with ahort time acale and high 



ent ropy, where intenaivc flux of free neutrona are ab — 
sorbed by seed elements successively to form the nu- 
clear reaction flow on extremely unstable nuclei in 
neutron-rich side. Recent progress in the studies of 
nuclear physics of unstable nuclei has made it possi- 
ble to simulate the r-process nucleosynthesis by the 
use of accumulated knowledge on nuclear masses and 
beta half-lives of several critical radioactive elements. 

The studies of r-process elements make another im- 
pact on the cosmic age problem, that is the age of 
the Universe to be known from cosmological constants 
and the age of the oldest globular cluster conflict with 
each other. A typical r-process element, thorium, has 
been detected recently in very metal-deficient stars, 
providing independent method to estimate the age of 
the Milky Way Galaxy ( ^eden et al. 19961) . Since 
thorium has half-life of 14 Gyr, the observed abun- 
dance relative to the other stable elements is used 
as a chronometer dating the age of the Galaxy. To 
study the origin of the r-process elements is thus im- 
portant and even critical in cosmology and astron- 
omy of Galactic chemical evolution as well as nuclear 
physics of unstable nuclei. Unfortunately, however, 
astrophysical site of the r-process nucleosynthesis has 
been poorly known, although several candidate sites 
are proposed and being investigated theoretically. 

Neutrino-driven wind, which is our object to study 
in this article, is thought to be one of the most promis- 
ing candidates for the r-process nucleosynthesis. It is 
generally believed that a neutron star is formed as 
the remnant of gravitational core collapse of Type 



II, iftr 



ui ie supeinuvau. — The liuL neuLiun sLai jusL 
boi' u releases mosL of its eneig}' as ueutiiuos dmiug 
Kelvin-Helmholtz cooling phase, and these neutrino 
drive matter outflow from the surface. This outflow is 
called neutrino-driven wind. Many theoretical stud- 



1994 ) . A benchmark study of numerical simulation by 



Wilson and his collaborators (Woosley et al. 1994) can 
successfully explain the solar system r-process abun- 



dances, but the others (Witti, Janka, & Takahashi 



1994, Takahashi, Witti, & Janka 1994) can not re- 



produce their result. Qian and Woosley (1996) tried 
to work out this discrepancy using approximate meth- 
ods to solve the spherically symmetric, steady state 
flow in the Newtonian framework. 

They could not find suitable condition for the 
r-process nucleosynthesis, and they suggested in a 
post-Newtonian calculation that general relativistic 
effects may improve thermodynamic condition for 
the r-process nucleosynthesis. Cardall and Fuller 
(1997) adopted similar approximate methods in gen- 
eral relativistic framework and obtained short dy- 
namic time scale of the expansion and large entropy, 
which is in reasonable agreement with the result in 
post-Newtonian approximation adopted by Qian and 
Woosley (1996). They did not remark quantita- 
tively, however, what kind of specific effect among 
several general relativistic effects is responsible for 
this change. 

Since the wind blows near the surface of the neu- 
tron star, it is needed to study expansion dynam- 
ics of neutrino-driven wind in general relativity. The 
first purpose of this paper is to quantitatively make 
clear the effects of general relativity by adopting fully 
general relativistic framework. Although we assume 
only spherical steady-state flow of the neutrino-driven 
wind, we do not adopt approximate methods as in 
several previous studies. We try to extract as gen- 
eral properties as possible of the wind in manners in- 
dependent of supernova models so that they are to 
be compared with expansion of different object like 



accretion disk of binary n eutron star merger ( 3ym 



balisty fc Schramm 1982) or sub-criti cal small mass 
neutron star ( ^umiyoshi et al. 1998| ), which is in- 
duced by intense neutrino burst. The second purpose 
is to look for suitable conditions for the r-process. 
There are key quantities in order to explain the so- 
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lar system r-proccss abundances. They arc the mass Wasserman 1986, Qian & Woosley 1996, Cardall & 



outflow rate, A/, the dynamic time scale of the expan- 



sion, Tdyn, the entropy, S, and the electron fraction. 
Ye. The third purpose of this paper is to make clear 
how these thermodynamic and hydrodynamic quan- 
tities affect the r-process nucleosynthesis by carrying 
out the nucleosynthesis calculation numerically. 

In the next section we explain our theoretical mod- 
els of neutrino-driven wind. We introduce basic equa- 
tions to describe the dynamics of the wind in the 
Schwarzschild geometry. Boundary conditions and 
adopted parameters for solving these equations are 
presented in this section. Numerical results are shown 
in section 3, where the effects of general relativity are 
studied in detail. We also investigate the dependence 
of the key physical quantities like Tdyn and S on the 
neutron star mass, radius, and neutrino Imninosity in 
order to look for the conditions of the neutrino-driven 
wind which is suitable for the r-process nucleosynthe- 
sis. Applying the result obtained in section 3, we 
carry out the nucleosynthesis calculation in section 4. 
The purpose of this section is to confirm quantita- 
tively that the r-process elements are produced suc- 
cessfully in the wind having very short dynamic time 
scale with relatively low entropy. We finally sum- 
marize the results of this paper and present further 
discussions and outlook in section 5. 

2. Models of neutrino-driven winds 
2.1. Basic equations 

Type II or lb supernova explosion is one of the com- 
plex hydrodynamic process which needs careful theo- 
retical studies of the convection associated with shock 
propagation. The time of our interest, however, is the 
later phase after the core bounce, at which the shock 
has already passed away to reach a radius about 10000 
km and continuous mass outflow is installed from the 
surface of the neutron star. Recent three dimensional 
numerical simulation (Hillebrandt 1998) has indicated 
that the convection near the shock front does not grow 
as deep as that shown in two-dimensional numerical 
simulation and the hydrodynamic conditions behind 
the shock are more likely similar to those obtained 
in one-dimensional numerical simulation. Since Wil- 
son's numerical simulation of SN1987A in Woosley et 
al. (1994) has shown that the neutrino-driven wind is 
adequately described by a steady state flow, we here 
adopt spherically symmetric and steady state wind, 
following the previous studies (Duncan, Shapiro, & 



Fifller 1997). According to his numerical simulation. 



the neutrino luminosity L^, changes slowly from about 
Iq52 gj.ggy's below 10^^ ergs/s during ~ 10 s of the 
Kelvin-Helmholtz cooling phase of the neutron star. 
The properties of the protoneutron star, i.e. the mass 
M and radius R, also evolve slowly. We therefore 
take these quantities L^,, M, and R as input param- 
eters in order to describe more rapid evolution of the 
neutrino-driven wind. 

The basic equations to describe the spherically 
symmetric and steady state winds in Schwarzschild 
geometry are given by (Shapiro & Teukolsky 1983) 



M — Anr pi,u, 



du 

dr 



dP 



Plot 



P dr 



2M 



M 



de 
dr 



Pi 



(1) 



(2) 



(3) 



where M is the mass outflow rate, r is the distance 
from the center of the neutron star, is the baryon 
mass density, u is the radial component of the four 
velocity, ptot = Pb + Pb£ is total energy density, s is 
the specific internal energy, P is the pressure, M is 
the mass of the neutron star, and q is the net heating 
rate due to neutrino interactions with matter. We 
use the conventional units that the plank constant h, 
the speed of light c, the Boltzmann constant k, and 
gravitational constant G, are taken to be unity. Since 
the neutrino-driven wind blows from the surface of 
the hot protoneutron star at high temperature T ~ 
5 MeV and also the physics of the wind is mostly 
determined at T ^ 0.5 MeV dQian fc Woosley 1996| ), 
the equations of state are approximately written as 



P = 



180 ruN 
llTT^ 3 T 

60 ph 2 rrLN 



(4) 
(5) 



where T is the temperature of the system, and m^v 
is the nucleon rest mass. We have assumed that the 
material in the wind consists of photons, relativistic 
electrons and positrons, and non-relativistic free nu- 
cleons. 

The heating rate q in Eq. ^ through the in- 
teractions between neutrinos and material takes the 
key to understand the dynamics of the neutrino- 
driven wind. Following Bethe (1993) and Qian and 
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Woosley (1994), we take account of the following 
five neutrino processes; neutrino and antineutrino ab- 
sorption by free nucleons, neutrino and antineutrino 
scattering by electrons and positrons, and neutrino- 
antineutrino annihilation into electron-positron pair 
as the heating processes, and electron and positron 
capture by free nucleons, and electron-positron anni- 
hilation into neutrino-antineutrino pair as the cool- 
ing processes. We assume that neutrinos are emitted 
isotropically from the surface of the neutron star at 
the radius R, which proves to be a good approxima- 
tion iii_£eccn^jiurM neutrino trans- 
fer ( Yamada, Janka, fc Suzuki 199^ ). In this paper. 



therefore, we make an assumption that the neutri- 
nosphere radius is equal to the protoneutron star ra- 
dius Ri, — R. Since the neutrino trajectory is bent 
in the Schwarzschild geometry, the material in the 
wind sees neutrinos within the solid angle subtended 
by the neutrinosphere which is greater than the solid 
angle in the Newtonian geometry at the same coordi- 
nate radius. The bending effect of the neutrino tra- 
jectory increases the heating rate compared to New- 
tonian case. We have to take account of the redshift 
effect on the neutrino energy, too, which tends to de- 
crease the heating rate. 

The important heating rate is due to the neutrino 
and antineutrino absorption by free nucleons 



and it is given by 
gi « 9.65A^a[(1 



■ p + e 



(6) 
(7) 



(8) 



where the first and second terms in the parenthesis 
are for the processes ^ and (Q), respectively, £i is the 
energy m MeV defined by Si = ^< Ef > / < Ei >, 
and < E^ > denotes the nth energy moment of the 
neutrino {i = Ve) and antineutrino (i = i>e) energy 
distribution, Na is the Avogadro number, is the 
electron fraction, i^^si is the individual neutrino or 
antineutrino luminosity in units of 10^^ ergs/s, and 
R,jQ is the neutrinosphere radius in units of 10^ cm. In 
this equation, l — gi{T) is the geometrical factor which 
represents the effect of bending neutrino trajectory, 
and gi{r) is given by 



where the function {1-2M /r) / {1-2M / R^,) arises due 
to the Schwarzschild geometry, and unity should be 
substituted for this factor in the Newtonian geometry. 
We also define the redshift factor 



$(r) = 



1 - 2M/R^ 
1 - 2M/r ' 



(10) 



in the Schwarzschild geometry, which is unity in the 
Newtonian geometry. We will discuss the effects of 
these general relativistic correction factors in the next 
section. 

The second heating rate due to neutrino and an- 
tineutrino scattering by electrons and positrons plays 
equally important role. Neutrinos of all flavors can 
contribute to the scattering, and the heating rate is 
given by 



93 ~ 2.nNA 



PS 



X l^|lM$(r)5MeV s-ig-i. 



(11) 



where =< Ef > / < Ei > in MeV (i = 
Ve, Vf., and v^), and we have assumed the same 
contribution from i^^, i^^, Vt, and Pt fluxes. We take 

~ l-14e| from the numerical studies by Qian and 
Woosley (1996). 

The third heating rate due to neutrino-antineutrino 
pair annihilation into electron-positron pair is given 

by 



15 



12.0Na 



6 



X ( Li.,,5i-^Pe,5i(e!.e + epJ + 



PsRte 



$(r)^MeV s-ig~\ 



(12) 



where .92 (^) is given by 

g2{r) = (1 - gi{r)f{9i{rf + 4gi(r) + 5). (13) 

The cooling rates which we included in the present 
calculations are for the inverse reactions of the two 
heating processes considered in Eqs. (||) and (p^. The 
first cooling rate due to electron and positron captures 
by free nucleons, which are the inverse reactions of (^ 
and (0), is given by 



92 « 2.27NATl,ev^eY s'^g" 



(14) 
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The second cooling rate due to electron-positron pair 


pernovae ( 


Woosley et al. 1994 




Witti, Janka, & Taka- 


am ihilation into ncutrino-antincutrino pair of all fla- 


hashi 1994, Takahashi, Witti, & Janka 1994 


), the 



vors, which is the inverse reaction of Eq. (^_2|), is given 

by 

q4 « 0.144iVA^^MeV s'^g-^ (15) 
Ps 

Combining the above five heating and cooling rates, 
we obtain the total net heating rate q 



g = 91 - 92 + 93 - 94 + 95- 



(16) 



As we will discuss in the next section, the first three 
heating and cooling rates gi , (72 , and dominate over 
the other two contributions from q^ and (75. 

2.2. Boundary conditions and input parame- 
ters 

We assume that the wind starts from the surface 
of the protoneutron star at the radius = i? and the 
temperature T^. Near the neutrinosphcre and the neu- 
tron star surface, both heating (mostly qi) and cool- 
ing (mostly 92) processes almost balance with each 
other due to very efficient neutrino interactions with 
material. The system is thus in kinetic equilibrium 
(Barrows and Mazurek 1982) at high temperature and 
high density. The inner boundary temperature Ti is 
determined so that the net heating rate q becomes 
zero at this radius. We have confirmed quantitatively 
that a small change in Ti does not influence the cal- 
culated thermodynamic and hydrodynamic quantities 
of the neutrino-driven wind very much. We give the 
density p{ri) — 10^*^ g/cm'^ at the inner boundary, 
which is taken from the result of Wilson's numerical 
simulation in Woosley et al. (1994). 

The luminosity of each type of neutrino Li (i — 
i/e, Ve, J^/i, i^fi, i^rj i^t) is similar to one another and 
changes from about 10^^ to lO'^" ergs/s very slowly 



during ^ 10 s ( Woosley et al. 1994 ). We therefore 
take a common neutrino luminosity as a constant 
input parameter. In the heating and cooling rates, 
however, we use the values of neutrino energies Cj^^ = 
12 MeV, e^- = 22 MeV, and = = 34 MeV for the 
other flavors at = i? as in Qian & Woosley(1996). 
We take the neutron star mass as a constant input 
parameter ranging 1.2M© < M < 2.OM0, too. 

The mass outflow rate M determines how much 
iorial ie ojoctod by tho neutrino driven wind. — 



Eq a.([^) (H), M ie taken to bo a conetant value to bo 
determined by the following outer boundary condi- 
tion. In any delayed explosion models of Type II su- 



shock wave moves away at the radius around 10000 
km above the neutron star surface at times Is <J t 
after the core bounce. As we stated in the previous 
subsection, the neutrino-driven wind is described by a 
steady state flow fairly well between the neutron star 
surface and the shock. From this observation, a typi- 
cal temperature at the location of the shock wave can 
be used as an outer boundary condition. We impose 
the boundary condition only for subsonic solutions by 
choosing the value of M < Merit so that T = O.lMeV 
at r ~10, 000km, where Merit is the critical value for 
supersonic solution. Given p(ri), Eq.(|^) determines 
also the initial velocity at r = for each M. 

We here explore the effects of the assumed bound- 
ary condition and the mass outflow rate M on the 
results of calculated quantities of the neutrino-driven 
winds. We show in Figs. 0(a) and |](b) the fluid veloc- 
ity and the temperature as a function of raduis from 
the center of neutron star for various M, where neu- 
tron star mass M = 1.4 Mq and neutrino luminosity 
Lv^ = 10^^ ergs/s are used. Figures ||(a) and ||(b) 
are the same as those in Figs. |^(a) and Qb) for M 
= 2.OM0 and L^^ = 10^^ ergs/s. Varied M's are tab- 
ulated in Table 1 with the calculated entropies and 
dynamic timescales. These figures indicate that both 
velocity and temperature profiles are very sensitive to 
the adopted M corresponding to different boundary 
conditions at r = 10000 km. However, the entropies 
are more or less similar to one another, while exhibit- 
ing very different dynamic timescales. 

Although finding an appropriate boundary condi- 
tion is not easy, it is one of preferable manners to 
match the condition obtained in numerical simula- 
tions of the supernova explosion. We studied one of 
the successful simulations of 2OM0 supernova explo- 
sion assuming M = I.4M0 (Wilson 1998). Exten- 
sive studies of the r- process (Woosley et al. 1994) 
are based on his supernova model. Careful obser- 
vation tells us that, although the neutrino luminos- 
ity for each flavor changes from 5 x 10^^ ergs/s to 
lO^^ergs/s, the temperature lowers progressively to 
0.1 MeV around r = 10,000 km where the shock front 
almost stays during 10 s after the core bounce at 
times which we are most interested in. It is to be 



noted that for successful r-process (Woosley et al 



1994) the temperature has to decrease gradually down 
to around 0.1 MeV at the external region. This will 
be discussed in later sections. As displayed in Figs. 
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||(a) and |i|(b), our calculation denoted by "3" meets 
with this imposed boundary condition. Although it 
may not be necessarily clear, we can adopt the same 



boi ndary condition for different neutron star masses 



which we study in this article, expecting that the 
physics continuously changes and also aiming at com- 
paring the results with one another which arise from 
the same boundary condition. Even in the case of 
massive neutron star having M = 2.0Mq, as displayed 
in Figs. |^(a) and |^(b), we can still find a solution de- 
noted by " 1" which satisfies the same outer boundary 
condition. Although we fortunately found a solution 
with reasonable value of M, careful studies of the nu- 
merical simulation in the case of massive neutron stars 
are highly desirable in order to find better boundary 
condition. 

Let us discuss how our adopted outer boundary 
condition is not unresonable. We are interested in 
the times Is < t when the neutrino-driven wind be- 
comes quasi steady state flow between the neutri- 
nosphere and the shock front. Intense flux of neu- 
trinos from the hot proto-neutron star have already 
interacted efficiently with radiation and relativistic 
electron-positron pairs at high temperature. Thus we 
have T ^ T^, where T and T^, are respectively the 
photon and neutrino temperatures. In this stage, the 
gain radius Rg ( Bcthe fc Wilson 1985| ) at which the 
neutrino heating and cooling balance with each other 
is very close to the neutrinosphere. Since we make an 
approximation that the neutrinosphere and the neu- 
tron star surface is close enough, we here assume that 
the gain radius is also the same, i.e. Rg — Rj^ — R. 
On these conditions we can estimate the mass out- 
flow rate M by considering the energy deposition to 
the gas from the main processes of neutrino capture 
on nucleons (||) and (^. 

Following the discussion by Woosley et al. (1994), 
the rate of energy deposition in the gas above the 
neutrinosphere is given by 



E = {L^^ + LpJ X Ti,, 



(17) 



where t^, is the optical depth for the processes (||) and 
(j^) and is given in terms of the opacity Ki, and the 
pressure scale height Lp by. 



Ky{Rg)p^{Rg)Lp{Rg) 



omeR^ 



3.5MeV 



1AM, 







M 



(18) 



Note that Rg = R and = T^. In order to ob- 
tain this expression, we have already used an ap- 
proximate opacity ( Bethe 199C| , [Woosley fc Weaver 
"199^ ) « 6.9 X 10^^«(r^/3.5MeV)^cm^g-^ and the 



pressure scale height in radiation dominated domain 
which is written as, 



(aT^)/{GMp^/R') 
74km 



(ilv) 



(19) 



where the subscripts on i?7 and pb,? indicate cgs mul- 
tipliers in units of 10^. The energy deposition Eq.(pT|) 
is mostly used for lifting the matter out of the grav- 
itational well of the neutron star. Thus, inserting 
Eq. (|l8|) into Eq. ( pT[ ) and using the relation L,y^ = L^^ 
= {7/4)nR'^aTI^, the mass outflow rate M is approx- 
imately given by 



M 



E/{GM/R) 



0.092 



ylO^^ergs s" 



5/2 



I AM, 



M 



Mr.,S- 



(20) 



Our mass outflow rate M obtained from the imposed 
boundary condition of a temperature 0.1 MeV at 
10,000 km is in reasonable agreement with the es- 
timate using this Eq. ( pO| ) within a factor of flve for 
lO^^ergs/s < (L^^ -1- L^-J < 1052gj.gs/g^ 

2.3. Characteristics of the neutrino-driven wind 

When the material of the wind is on the surface of 
the neutron star and neutrinosphere, thermodynamic 
quantities still reflect the effects of neutralization and 
the electron fraction Yf. remains as low as ^ 0.1. Once 
the wind leaves surface after the core bounce, electron 
number density decreases abruptly and the chemical 
equilibrium among leptons is determined by the bal- 
ance between the two processes (^) and (|^) due to 
intense neutrino fluxes, shifting to ~ 0.5. Interest- 
ing phase starts when the temperature falls to ^ 10^° 
K, for our purpose of studying the physical condition 
of the neutrino-driven wind that is suitable for the 
r-process nucleosynthesis. At this temperature the 
material is still in the NSE, and the baryon numbers 
are carried by only free protons and neutrons. The 
neutron-to-proton number abundance ratio is deter- 
mined by Yg for charge neutrality. 
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Electron antineutrino has a harder spectrum than 
electron neutrino, as evident from their energy mo- 
ments ty^ = 12 MeV < ep^ = 22 MeV. Thus, the 
material is slightly shifted to neutron-rich. Assum- 
ing weak equilibrium, this situation is approximately 
described by 



1 



25+ 1.2(52/ep 



(21) 



where Ajy_,„ and Ap^p are the reaction rates for the pro- 
cesses (H) and (0), resp ectively, and b is the n eutron- 
proton mass difference ( Qian fc Woosley 1996 ). In our 
parameter set of the neutron star mass M = 1.4Mq 
and radius i? = 10 km, for example, Kg varies from 
Y^(r = R) = 0.43 to y;(r = 10000km) = 0.46 very 
slowly due to the redshift factor ( [lO| ) because of e oc 
As this change is small and the calculated result of 
hydrodynamic quantities are insensitive to Yg, we set 
Ye = 0.5 for numerical simplicity. 

One of the most important hydrodynamic quan- 
tity, that characterizes the expansion dynamics of the 
neutrino-driven wind, is the dynamic time scale Tdyn 
which is the duration time of the a-process. When 
the temperature falls below 10^'^ K, the NSE favors 
a composition of alpha-particles and neutrons. As 
the temperature drops further below about 5 x 10^ 
K (T « 0.5 MeV), the system falls out of the NSE 
and the a-process starts accumulating some amount 
of seed elements until the charged particle reactions 
freeze out at T « 0.5/e MeV « 0.2 MeV. Introduc- 
ing a time variable of the wind moving away from the 
distance to outer distance r/ 

'^^ dr 

(22) 



and setting = r{T = 0.5MeV) and r/ = r(T = 
0.5/e MeV), we can define the dynamic time scale 
Tdyn by 



T=0.5/o McV 



dr 



Tdyn = I ■ (23) 

/t=0.5McV " 

The second important hydrodynamic quantity, that 
affects strongly the r-process nucleosynthesis which 
occurs at later times when the temperature cools be- 
low 0.2 MeV, is the entropy per baryon, defined by 



uT 



dr. 



(24) 



where q is the total net heating rate (p^). As 5 oc 
/ Pb assuming the radiation dominance, high en- 
tropy and high temperature characterizes a system 
with many photons and low baryon number density. 
Since high entropy favors also a large fraction of free 
nucleons in the limit of the NSE, it is expected to be 
an ideal condition for making high neutron-to-seed 
abundance ratio. Therefore, the high entropy at the 
beginning of the a-process is presumed to be desirable 
for successful r-process. 

3. Numerical results 

3.1. Effects of relativistic gravity to entropy 

The purpose of this section is to discuss both sim- 
ilarities and differences of the neutrino-driven wind 
between the relativistic treatment and the Newtonian 
treatment. In Fig. we show typical numerical re- 
sults of radial velocity u, temperature T, and baryon 
mass density pb of the wind for the neutron star mass 
M = 1.4M0, radius R = 10km, and the neutrino lu- 
minosity Ly = 10^"'^ ergs/s. The radial dependence 
of these quantities is displayed by solid and dashed 
curves for Schwarzschild and Newtonian cases, respec- 
tively, in this figure. Using these results and Eq. (p^) , 
we can calculate S in each ejecta. Figure ^ shows the 
calculated profile of the entropy S for the two cases. 
Although both entropies describe rapid increase just 
above the surface of the neutron star 10 km < r < 15 
km, the asymptotic value in general relativistic wind 
is nearly 40 % larger than that in Newtonian wind. 

The similar behavior of rapid increase in both 
winds is due to efficient neutrino heating near the 
surface of the neutron star. We show the radial de- 
pendence of the heating and cooling rates by neutri- 
nos in Figs. ||(a)-(c). Figure |^(a) shows the total 
net heating rate defined by Eq. (p6[), and Figs. |^(b) 
and (c) display the decompositions into contribution 
from each heating(solid) or cooling(dashed) rate in 
Schwarzschild and Newtonian cases, respectively. The 
common characteristic in both cases is that net heat- 
ing rate q has a peak around r « 12 km, which makes 
a rapid increase in S near the surface of the neutron 
star for the following reason. The integrand of the 
entropy S in Eq. (|2j) consists of the heating rate and 
the inverse of fluid velocity times temperature. The 
fluid velocity increases more rapidly than the slower 
decrease in the temperature, as shown in Fig. |^, after 
the wind lifts off the surface of the neutron star. 

Let us carefully discuss the reason why the general 
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relativistic wind results in 40 % larger entropy than 
the Newtonian wind in the asymptotic region. This 
fact has been suggested in the previous papers of Qian 
k Woosley (1996) and Cardall & Fuller (1997). Un- 
fortunately, however, the reason of this difference was 
not clearly appreciated to the specific effect quantita- 
tively among several possible sources. 

We first consider the redshift effect and the bend- 
ing effect of the neutrino trajectory. The redshift 
effect plays a role in decreasing the mean neutrino 
energy e,y ejected from the neutrinosphere, and in 
practice €i, is proportional to the redshift factor <i>(r) 
which is defined by Eq. (^0|) . Since neutrino luminos- 
ity is proportional to <i>^ and the heating rate (71,93, 
and <75 depend on these quantities in different man- 
ners, each heating rate has different ^-dependence as 
qi oc L^el oc oc L^e^, oc and 55 cx i^e^ oc 

$9, as shown in Eqs. (|),(0), and (|l|). Coohng 
rates (72 and (74 do not depend on <I>(r). The bending 
effect of the neutrino trajectory is included in the ge- 
ometrical factors gi{r) and 52 (^) in these equations. 
Although numerical calculations were carried out by 
including all five heating and cooling processes, as (71, 
q2, and (73 predominate the total net heating rate q, 
we here discuss only these three processes in the fol- 
lowing discussions for simplicity. 

In Newtonian analysis, the redshift factor <i>(r) is 
unity and the geometrical factor is given by 



9iN{r) 



This geometrical factor gi{r) and the redshift factor 
appear in the form of (1 — gi(r))$(r)'" in the heating 
rate qi {m — 5) and (73 (m = 6). As for the first fac- 
tor (1 — (7i(r)), the following inequality relation holds 
between the Schwarzschild and Newtonian cases, for 
Rv < r; 

(l-(7i(r)) > (l-gi^(r)). 

However, $(r) is a monotonously decreasing func- 
tion of r, the combined factor (1 — (7i(r))$(r)'"/(l ~ 
giNif)) increases from unity and has a local max- 
imum around r ~ Ri, + 0.2 km. Its departure from 
unity is at most 3 % . Beyond this radius the function 
starts decreasing rapidly because of the redshift effect 
$(r)'", and it becomes as low as ~ 0.6 at r ~ 30 km. 
In this region, the net heating rate in the relativistic 
wind is smaller than that in the Newtonian wind if the 
temperature and density are the same. However, the 



difference in this region does not influence the dynam- 
ics of the wind very much. It is almost determined in 
the inner region i?^ < < 15 km where one finds ef- 
ficient neutrino heating and small difference between 
(l-(7i(r))$™(r) and (l-g.^ir)). 

By performing general relativistic calculation and 
neglecting these two relativistic effects, i.e. the red- 
shift effect and the bending effect of the neutrino tra- 
jectory, we find that it produces only a small change 
in entropy by AS* ^ 3. Thus it does not seem to be 
the major source of the increase in the entropy. 

Let us consider another source of general relativis- 
tic effects which are included in the solution of a set 
of the basic equations (0)-(||)- Since the entropy de- 
pends on three hydrodynamic quantities q{r), u{r), 
and T(r) (see Eq.(p4[)), we should discuss each quan- 
tity. The neutrino- heating rate, q{r), depends on the 
temperature T(r) and density phir) in addition to 
the redshift factor and the geometrical factor of the 
bending neutrino trajectory. Therefore, we study first 
the detailed behavior of T{r), u{r), and Phif), and 
then try to look for the reason why the general rela- 
tivistic effects increase the entropy. We assume that 
the pressure and internal energy per baryon are ap- 
proximately described by the radiation and relativis- 
tic electrons and positrons in order to make clear the 
following discussions. This is a good approximation 
for the neutrino-driven wind. The equations of state 
are given by 

P^'-^T', (25) 

_ llTT^T^ 

By using another approximation 



(26) 



4F 



(27) 



which is satisfied in the region of interest, we find 

]_dT _ 1 ph + P 

T dr 1 + u2 _ 2M 4p 

r 

M 2m2 45 upb 
72" + ~ ~ 117^2 T4 ' 



(28) 



in Schwarzschild case. The basic equations of the 
spherically symmetric and steady state wind in New- 
tonian case are given by 



M — Anr p^v, 



(29) 
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dv 



1 dP M 
72"' 



dr Ply dr 1 

de P dpb 
^ dr dr 



(30) 



(31) 



where v is the fluid velocity. The equations of state 
are given by Eqs. (jj) and (||) the same as in Schwarzschild 
case. Repeating the same mathematical technique in 
Eqs. (^-(|l|) instead of Eqs. (|5)-(|) and taking the 
same approximations as ( |25|)-(|2^ ), we find the equa- 
tion corresponding to Eq. (^8|), in Newtonian case, 
as 



1 dT 



Pb_ 

AP 



M 



r 



45 vph 
Il7r2 7^' 



(32) 



Note that the logarithmic derivative of the temper- 
ature, dlnT/dr = T~^dT/dr, has always a nega- 
tive value, and the temperature is a monotonously 
decreasing function of r. There are two differences 
between Eqs. (28) and (|3^). The first prefactor 
1/(1 -f - 2M/r) in the r.h.s. of Eq. (|2|) is larger 
than unity. This causes more rapid decrease of T(r) 
in relativistic case than in Newtonian case at small 
radii within r 20 km, as shown in Fig. ^, where 
our approximations are satisfied. The second prefac- 
tor (pb + P)/4P in the r.h.s. of Eq. (|2|) is larger 
than the prefactor Ph/^P in the r.h.s. of Eq. ([3^), 
i.e. (pb + P)/4:P > pb/AP, which also makes the 
difference caused by the first prefactor even larger. 

Applying the similar mathematical transformations 
to the velocity, we obtain the following approxima- 
tions, 



1 du 
u dr 



(pb + 4P) M 2 



Pb 



i + u^-m AP 

r 

in Schwarzschild case, and 

Idv ^ 3pb M 2 
V dr AP 3r 



Pb ^ 

AvP^ 



(34) 



in Newtonian case. In these two equations, the first 
leading term in the r.h.s. makes the major contribu- 
tion. Since exactly the same prefactors 1/(1 + — 
2M/r) and [ph + AP)/AP appear in Schwarzschild 
case, the same logic as in the logarithmic derivative of 
the temperature is applied to the velocity. Note, how- 
ever, that slightly different initial velocities at the sur- 
face of the neutron star make this difference unclear 
in Fig. H The relativistic Schwarzschild wind starts 
from u(10 km) « 8.1 x 10"' cm/s, while the Newtonian 
wind starts from v{10 km) w 2.0 x 10^ cm/s. Both 



winds reach almost the same velocity around r ~ 20 
km or beyond. 

The baryon number conservation leads to the log- 
arithmic derivative of the baryon density 



1 dpb 
Pb dr 



1 du 
u dr 



(35) 



where u is the radial component of the four- velocity in 
Schwarzschild case. The fluid velocity v should read 
for u in Newtonian case. Inserting Eq. ( |3^ ) or Eq. (34) 
to the first leading term of the r.h.s. of this equation, 
we can predict the behavior of pb as a function of r 
in both Schwarzschild and Newtonian cases as shown 
in Fig. I 

Incorporating these findings concerning T(r), and 
u(r) into the definition of entropy Eq. ([24|), we can 
now discuss why the relativistic Schwarzschild wind 
makes larger entropy than the Newtonian wind. We 
have already discussed previously in the second para- 
graph of this section that the fluid velocity increases 
more rapidly in Schwarzschild case. Since integrand 
of the entropy S is inversely proportional to fluid ve- 
locity times temperature, this fact enlarges the dif- 
ference due to q at smaller radii (see Fig. ^(a)). In 
addition, as we found, temperature in Schwarzschild 
geometry is smaller than the temperature in Newto- 
nian geometry. For these reasons, the entropy in the 
relativistic Schwarzschild wind becomes larger than 
the entropy of the Newtonian wind. 

Let us confirm the present results quantitatively 
in a different manner. The entropy per baryon for 
relativistic particles with zero chemical potential is 
given by 

S.H^4^. (36) 

45 Pb/rriN 

Here, we take a common temperature T — 0.5 MeV 
to each other in Schwarzschild and Newtonian cases. 
This is the typical temperature at the beginning of 
the a-process, and both electrons and positrons are 
still relativistic at this temperature. We read off the 
radii at which the temperature becomes 0.5MeV in 
Fig. H They are 43km and 55km in Schwarzschild 
and Newtonian cases, respectively. We can again read 
off the baryon mass densities at these radii in this fig- 
ure, that are pb — 5.5 x lO'^ g/cm'^ at r = 43 km 
in relativistic Schwarzschild wind and pb = 7.8 x 10^ 
g/crn^ at r = 55 km in Newtonian wind. Taking the 
inverse ratio of these pb values with approximate re- 
lation (am , we find that the entropy in Schwarzschild 
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case is 40 % larger than that in Newtonian case. This 
is quantitatively in good agreement with the result of 
numerical calculation shown in Fig. ^. 

Let us shortly remark on the dynamic time scale 
Tdyn- Although higher entropy is favorable for making 
enough neutrons in the neutrino-driven wind, shorter 
dynamic time scale also is in favor of the r-process. 
This is because the neutron-to-seed abundance ratio, 
which is one of the critical parameters for success- 
ful r-process, becomes larger in the wind with shorter 
Tdyn, which is to bc discussed in the next section. 
It is therefore worth while discussing the general rel- 
ativistic effect on Tdyn here. The argument is very 
transparent by using Eqs. and (|3^) and Fig. |^. 
Since the dynamic time scale Tdyn is defined as the du- 
ration of a-process in which the temperature of the 
wind cools from T = 0.5 MeV to T = 0.5/e ss 0.2 
MeV, faster cooling is likely to result in shorter Tdyn- 
Let us demonstrate it numerically. For the reasons 
discussed below two Eqs. (p^ ) and js^), the relativis- 
tic fiuid describes more rapid decrease in temperature 
than the Newtonian fluid as a function of distance r. 
In fact, the distance corresponding to T = 0.5 — 0.2 
MeV are r = 43 — 192 km in Schwarzschild case, and 
r = 55 — 250 km in Newtonian case. Figure ^ tells 
us that both fluids have almost the same velocities 
at these distances, which gives shorter Tdyn for the 
Schwarzschild case than the Newtonian case. The cal- 
culated dynamic time scales are Tdyn = 0.164 s for the 
former and Tdyn = 0.213 s for the latter. 

Before closing this subsection, let us briefly discuss 
how the system makes a complicated response to the 
change in T(r), u{r) and Phif)- When the temper- 
ature decreases rapidly at 10 km < r < 20 km, the 
major cooling process of the e'^e~ capture by free 
nucleons, q2, is suppressed because this cooling rate 
has rather strong temperature dependence, q2 cx T^. 
In Schwarzschild geometry this suppression partially 
offsets the decrease in qi due to the neutrino redshift 
effect, though being independent of temperature of 
the wind. Another heating source q^ due to neutrino- 
electron scattering also plays a role in the change of 
entropy. Since q^ depends on the baryon density as 
well as temperature q^ oc T'^/p\^, if the system has 
a correlated response to decrease pb strongly with 
decreasing temperature, then this might eventually 
work for the partial increase in entropy. However, in 
reality, actual response arises from more complicated 
machanism because g^'s should depend on the solu- 
tion of dynamic equations (l)-(3) self-consistcntly on 



adopted proper boundary conditions and input pa- 
rameters through the relation qi oc Lj/eJ, q2 oc T^, 
qs oc T^/phL^e^, q^ oc T^/pb, and 55 oc p^^L^e^. 
The neutrino-driven wind is a highly non-linear sys- 
tem. 

3.2. Parameter dependence 

Most of the previous studies of the neutrino-driven 
wind have been concentrated on SN1987A, and the 
parameter set in the theoretical calculations was al- 
most exclusive. We here expand our parameter region 
of the neutron star mass M, radius R, and neutrino 
luminosity L^, and investigate widely the dependence 
of key quantities, Tdyn and S, on these three param- 
eters. Since the neutron star mass M and radius R 
are mostly contained through the form M/R in the 
basic equations of the system, we only look at the 
dependence on M and L^. 

Figures ||(a) and ||(b) show the calculated Tdyn 
and S at the beginning of the a-process at T = 0.5 
MeV for various neutron star masses I.2M0 < M < 
2.OM0. Closed circles, connected by thick solid line, 
and open triangles, connected by thin solid line, are 
those for the Schwarzschild and Newtonian cases. In 
Fig. H(a) , we plot also two broken lines in Newtonian 
case from the paper (Qian and Woosley 1996) which 
adopted 



Tdyn(QW) 



O.SMcV 



(37) 



in two limits of the radiation dominance (upper) and 
the dominance of non-relativistic nucleon (lower). In 
either limit, this Tdyn(QW) is an increasing function of 
the neutron star mass and this feature is in reasonable 
agreement with our exact solution Eq.(^). However, 
absolute value of (^)is about half that of the exact 
solution in the Newtonian case. 

Remarkable difference between Schwarzschild and 
Newtonian cases is an opposite response of Tdyn to 
the neutron star mass (Fig. ^(a)). General relativis- 
tic effects make the dynamic time scale even smaller 
with increasing neutron star mass. We have already 
discussed the reason why Tdyn in Schwarzschild case 
is smaller than that in Newtonian case by compar- 
ing Eqs. (|2^) and (|3^) from each other. We under- 
stand the decrease of Tdyn as a consequence from the 
fact that the general relativistic effects, which arise 
from the two prefactors in the r.h.s. of Eq.([28|), are 
enlarged by stronger gravitational force Mjr^ with 
larger M. Similar analysis on the role of the gravita- 
tional force is applied to the discussion of entropy and 
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Eqs. (11), (HI), and (H). Figure |(b) displays that 
the entropy per baryon in Schwarzschild case makes 
stronger mass dependence than in Newtonian case. 

It is to be noted again that the above features of 
the mass dependence are equivalent to those obtained 
by the change in the neutron star radius. Since the 
radius of protoncutron star shrinks with time in cool- 
ing process, it may work for increasing the entropy 
and decreasing the dynamic time scale. 

Figures ^(a) and 0(b) show the dependence of our 
calculated Tdyn and S on the neutrino luminosity 
ranging 10^'^ergs/s < < lO^^ergs/s. Differing from 
the mass dependence, both quantities are decreasing 
function of Li, as far as Li, < 10^^ ergs/s. This 
tendency, except for the absolute values, is in rea- 
sonable agreement with approximate estimates (Qian 
and Woosley 1996) shown by broken lines. This is be- 
cause larger luminosity makes the mass outflow rate 
M higher through more efficient neutrino heating, 
which causes bigger increase in the fluid velocity in 
addition to moderate increase in baryon density. Hav- 
ing these changes in hydrodynamic quantities with 
the definition of Tdyn, Eq.(|2^), and the definition of 
S, Eq.([2^), we understand that both quantities de- 
crease with increasing neutrino luminosity. 

However, if the luminosity becomes larger than 
]^g52 ergs/s, the temperature does not decrease as low 
as 0.1 MeV before the distance reaches 10000 km be- 
cause of the effect of too strong neutrino heating. The 
dynamic time scale Tdyn is of order ~ 10 s. In such a 
very slow expansion of the neutrino-driven wind, a- 
process goes on and leads to uninteresting r-process 
nucleosynthesis. 

To summarize this section, we find it difficult to ob- 
tain very large entropy 400 for reasonably short dy- 
namic time scale Tdyn < 0.1 s, as reported by Woosley 
et al. (1994), by changing the neutron star mass M 
and neutrino luminosity L^. However, there are still 
significant differences between our calculated result 
of Tdyn and S, which are shown by thick solid lines 
in Figs. ||(a)-^(b), and those of Qian and Woosley 
(1996), which arc shown by broken lines, in the mass 
dependence of the entropy and the opposite behavior 
in Tdyn. We will see in the subsequent sections that 
these differences are important to look for successful 
condition of the r-process. 



3.3. Implication in nucleosynthesis 

Having known the detailed behavior of dynamic 
time scale Tdyn and entropy per baryon S* as a func- 
tion of neutron star mass M, radius R, and neutrino 
luminosity L^, we are forced to discuss their impli- 
cation in the r-process nucleosynthesis. We have al- 
ready shown the calculated results of Tdyn and S for 
limited sets of two independent parameters M and 
Li, in Figs. ||(a)-0(b). We here expand the parameter 
space in order to include a number of (Af, Li^)-grids 
in their reasonable range 1.2Mq < M < 2.0Mq and 
lO^^ergs/s <L^< 10^^ ergs/s. 

Figure ^ displays the calculated results in the 
Tdyyi — S plane. Shown also are two zones for which 
the r-process nucleosynthesis might occur so that the 
second abundance peak around A — HQ and the third 
abundance peak around yl = 195 emerge from a the- 
oretical calculation as suggested by Hoffman et al. 
(1997). Their condition for the element with mass 
number A to be produced in an explosive r-process 
nucleosynthesis, for Yf. > {Z)/{A), is given by 

^^y \ 8xlO^((A)-2(Z)) /T,,„x)^/^ 
^ \ ln[(l - 2(Z)/A)/(1 - {A)IA)\ \sec)] ' 

(38) 

where (A) is mean mass number and (Z) is mean 
proton number of the seed nuclei at the end of the 
a-process. Following numerical survey of seed abun- 
dance of Hoffman et al. (1997), we choose (A) = 90 
and (Z) = 34 in Fig. ^. From this figure, we find 
that dynamic time scale as short as Tdyn ~ 6 ms with 
M =2.0Mq and = lO^^ergs/s is the best case 
among those studied in the present paper in order to 
produce the r-process elements, although the entropy 
S is rather small 140. 

Let us remark shortly on this useful equation. 
Equation (|3|) tells us that the r-process element with 
mass number A is efficiently produced from seed ele- 
ments with (A) and (Z) on a given physical condition 
Tdym S, and Ye at the onset of r-process nucleosynthe- 
sis at Tg w 2.5. In order to derive Eq. (|3^), Hoffman 
et al. (1997) assumed that the a + a + n -^^ Be -|- 7 
reaction is in equilibrium, because of its low Q-value, 
during the a-process at T w 0.5 — 0.2 MeV and that 
the ^Be -I- a — C -t- 7 reaction triggers burning 
of alpha-particles to accumulate seed elements. The 
NSE holds true if the nuclear interaction time scale for 
a + a + n — >^ Be -|- 7 is much shorter than the expan- 
sion time scale. We found in the present calculation 
that it is not always the case in neutrino-driven winds 
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with short dynamic time scale, for « 5 x 10^^ — 10^^ 
ergs/s, which is to be discussed more quantitatively in 
the next section. Keeping this in mind, we think that 
Eq. ( |3^ ) is still a useful formula in order to search for 
suitable physical condition for the r-process without 
performing numerical nucleosynthesis calculation. 

One might wonder if the dynamic time scale Tdyn ~ 
6 ms is too short for the wind to be heated by neu- 
trinos. Careful comparison between proper expan- 
sion time and specific collision time for the neutrino 
heating is needed in order to answer this question. 
Note that Tdyn was defined as the duration of the 
a-process so that the temperature of the expanding 
wind decreases from T = 0.5 MeV to 0.5/e « 0.2 
MeV, which correspond to outer atmosphere of the 
neutron star. These radii are r(T = 0.5 MeV) = 52 
km and r(T = 0.5/e MeV) = 101 km for the wind 
with {L^,M) = (1052gj.gs/s 2.0AfQ), and r{T = 
0.5 MeV) = 43 km and r{T = 0.5/e MeV) = 192 
km for the wind with {L^,M) = (lO^^^erg/s, I.4M0). 
We found in Figs. ||(a)-(c) that the neutrinos transfer 
their kinetic energy to the wind most effectively just 
above the neutron star surface at 10km < r < 20km. 
Therefore, as for the heating problem, one should re- 
fer the duration of time for the wind to reach the 
radius where temperature is T w 0.5 MeV rather 
than Tdyn- We can estimate this expansion time 
Theat by setting = i? = 10 km and rf — r{T ~ 
0.5MeV) in Eq. (|2|): Thcat = 0.017 s and 0.28 s for 
the winds with {L^,M) = (lO^^ergs/s, 2.OA/0) and 
(lO'^^ergs/s, 1.4Af0), respectively. We note, for com- 
pleteness, r(T = 0.5MeV) = 52 km or 43 km for each 
case. 

These proper expansion time scales, Thcat, are to 
be .cnmpfirpd with the snecifir cnllisinn time t„i fi 



thel neutrino-nucleus interartinns in order tn discuss 
the efficiency of the neutrino heating. The collision 



time Ti, is expressed (Qian et al. 1997) as 



0.201 X L 



u,51 



MeV J V 100km 



10- 



s(39) 



where -Li/,51 and e,y have already been defined in 
Sec. 2-1, and (cr^) is the averaged cross section over 
neutrino energy spectrum. As discussed above, neu- 
trino heating occurs most effectively at r « 12 km 
(see also Fig. ^(a)), and we set this value in Eq. (|39|). 
Since two neutrino processes and (0) make the 
biggest contribution to heating the wind and e^,^ — 12 



MeV and ep^ = 22 MeV, we set = {e^^ +epJ/2 ~ 15 
MeV. We take (cr^) = 10"''^ cm^. Incorporating these 
values into Eq. (|3^), we can obtain t^, value. Let us 
compare the specific collision time, t^, and the proper 
expansion time, Theatj with each other: 



0.0043s < Theat = 0.017 s, 

for (ii,,M) = (lO^^ergs/s,2.OM0), 



(40a) 

T^ = 0.043s < Theat = 0.28s, 

for {L^,M) ^ (lO^iergs/s,1.4M0). 

(40b) 

We can conclude that there is enough time for the 
expanding wind to be heated by neutrinos even with 
short dynamic time scale for the a-process, Tdyn ~ 6 
ms, which corresponds to the case (40a). 



Before closing this section, let us briefly discuss 
the effect of electron fraction Yg on the hydrodynamic 
condition of the neutrino-driven wind. Although we 
took Ye — 0.5 for simplicity in our numerical cal- 
culations, we should examine the sensitivity of the 
calculated result on quantitatively. Since we are 
interested in short dynamic time scale, let us inves- 
tigate the case with {L^,M) = (lO'^^ergs/s, 2.0Mq) 
which results in S* = 138.5 and Tdyn — 0.00618 s for 
Ye = 0.5. When we adopt Ye = 0.4, these quantities 
change slightly to S' = 141.5 and Tdyn = 0.00652 sec. 
These are very small changes less than 5%, and the 
situation in similar for the other sets of {L^, M). 

To summarize this section, we found that there is 
a parameter region in Fig. ^ which leads to desirable 
physical condition for the r-process nucleosynthesis. 



Sophisticated supernova simulation (Woosley et al 



1994) indicates that the neutrino luminosity from the 
protoneutron star decreases slowly from about 5x 10^^ 
to 10^^ ergs/s as the time passes by after the core 
bounce. Therefore, our favorable neutrino luminosity 
Li, — 10^^ ergs/s is possible in reality in relatively 
earlier epoch of supernova explosion at around 0.5 s 
to a few seconds after the core bounce. 

4. R-process nucleosynthesis calculation 

Our discussion on the r-process nucleosynthesis in 
the last section was based on Hoffman's criterion, 
Eq. (^8|), which is to be referred with caution for 
several assumptions and approximations adopted in 
its derivation. The purpose of this section is to con- 
firm quantitatively that the r-process occurs in the 
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neutrino-driven wind with short dynamic time scale, 
which we found in the present study. 

Given the flow trajectory characterized by u{t), 
Pb{t), and T{t) as discussed in the last section, our 
nucleosynthesis calculation starts from the time when 
the temperature is Tg = 9. Since this temperature is 
high enough for the system to be in the NSE, ini- 
tial nuclear composition consists of free neutrons and 
protons. We set Yf, — 0.4 in order to compare with 
Hoffman's criterion shown in Fig. |^. In our nucle- 
osynthesis calculation we used a fully implicit single 
network code for the a-process and r-process includ- 
ing over 3000 isotopes. We take the thermonuclear 
rea ction rates for all relevant nuclear processes and 



the ir inverse reactions as well as weak interactions 
from Thielemann (1995) for the isotopes Z < 46 and 
from Cowan et al. (1991) for the isotopes Z > 46. 
Previous r-process calculations had complexity that 
the seed abundance distribution at Tg — 2.5 was not 



this temperature are the baryon mass density ph = 
3.73 X lO''^ g/cm^, neutron mass fraction X„ = 0.159, 
mass fraction of alpha-particle = 0.693, average 
mass number of seed nuclei (A) — 94, and neutron- 
to-seed abundance ratio n/s — 99.8 for the set of 
hydrodynamic quantities Tdyn — 0.0062 s, S K 139, 
and Ye = 0.4. These values should be compared 
with those adopted in Woosley's calculation of tra- 
jectory 40, i,e, pb = 1.107 X 10" g/cm^, X„ = 0.176, 
Xa = 0.606, {A) = 95, n/s = 77, Tdyn w 0.305 s, 
S « 400, and Y^ = 0.3835, as in Table' 3 in Woosley 
et al. (1994). It is interesting to point out that our 
seed abun dance distribution in Fig. p| is very similar 
to th eirs ( Woosley et al. 1994 , Woosley fc Hoffman 
1992 ), as clearly shown by almost the same {A) w 95, 



full y shown in literature ([Woosley et al. 1994|, [Woosley 
& ] [loffman 19"92| , [Hoffman, Woosley, fc Qian 1997[ ), 
which makes the interpretation of the whole nucle- 
osynthesis process less transparent. This inconve- 
nience happened because it was numerically too heavy 
to run both a-process and r-process in a single net- 



although the other evolutionary parameters and ther- 
modynamic quantities are different from each other. 
The calculated final r-process abundance is displayed 
in Fig. Our wind model can produce the second 
{A « 135) and third {A « 195) r-process abundance 
peaks and rare earth elements between them as well. 

It is generally accepted that the r-process elements 
will be produced if there are plenty of free neutrons 
and if the neutron-to-seed abundance ratio is high 
enough to approximately sat isfy A « {A)+n/s (Hoff 



WO] k code for huge number of reaction couplings 



man, Woosley, fc Qian 1997 ) at the beginning of the 



among ~ 3000 isotopes. For this reason, one had 
to calculate the a-process first, using smaller network 
for light-to-intermediate mass elements, in order to 
provide seed abundance distribution at Tg = 2.5 (T w 
0.2MeV). Adopting such seed abundance distribution 
and following the evolution of material in the wind 
after T « 0.2 MeV, which is the onset temperature 
of the r-process, the r-process nucleosynthesis calcu- 
lation was extensively carried out by using another 
network code independent of the a-process. Our nu- 
cleosynthesis calculation is completely free from this 
complexity because we exploited single network code 
which is applied to a sequence of the whole processes 
of NSE - a-process - r-process. 

The calculated mass abundance distribution is shown 
in Figs. ^ and ^ for the neutrino-driven wind with 
{L^,M) = (10^2 ergs/s, 2. OA/©) that makes most fa- 
vorable condition for the r-process nucleosynthesis 
with the shortest Tdyn = 0.0062 s among those studied 
in the present paper (see Fig. ^ . Figure || displays 
the snapshot at the time when the temperature cooled 
to Tg = 2.5 (w 0.2MeV) at the end of the a-process. 
This shows seed abundance distribution at the onset 
of the r-process, too. Our calculated quantities at 



r-process, where A is the typical mass number of the 
r-process element. Therefore, the a-process should 
take the key to understand why our wind model re- 
sults in a similar r-process nucleosynthesis to the re- 
sult of Woosley's trajectory 40. 

The a burning starts when the temperature cools 
below T = 0.5 MeV. Since triple alpha reaction 
'''He(aa, 7)^^C is too slow at this temperature, alter- 
native nuclear reaction path to reach ^^C, ''He(an, 7) 
^Be(a,n)^^C, triggers explosive a-process to produce 
the seed elements. In rapidly expanding flow of 
neutrino-driven wind with short Tdyn, it is not a 
good approximation to assume that the first reaction 
"He(an,7)^Be is in the NSE. Rate equation is thus 
written as 



dYg 
dt 



plY^Y,A{aan Be) 

PbYcX^K^^c a C) 
[their inverse and 
other reaction rates), 



(41) 



where Ig, Yq,, Yn are the number fractions of ^Be, a, 
and neutron, and A(aan —^^ Be) and A(^Bea — C) 
are the thermonuclear reaction rate for each reaction 
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process as indicated. Details on A's are reported in 
Woosley and Hoffman (1992) and Wrean, Brune, and 
Kavanagh (1994). Let us talce the first term of the 
r.h.s. of Eq. ( p[ ) which is largest all terms in Eq. (|4l|). 
This is allowed in the following discussion of the time 
scale because the ^He(an, 7)^Be reaction is the slow- 
est among all charged particle reaction paths in all a- 
process reactions. We now define the typical nuclear 



rea:tion time scale of the a-process, regulated by 
the *Hc(Q!n, 7)''^Be reaction time scale tn, as 



(42) 



We show the ratio Tdyn/TN as a function of the baryon 
mass density pb at the beginning of the a-process 
when T = 0.5 MeV for various cases of the wind 
models with {L^ , M) in Fig. |l^. Note that the criti- 
cal line Tdyn/Ta = 1 is slightly shifted upwards be- 
cause of Tjv <^ Tq. This figure, with the help of 
Fig. H, clearly indicates that the favorable condi- 
tions for the r-process nucleosynthesis have inevitably 
shorter Tdyn ^ tjv and Tq,. Typical ratio is of order 
Tdyn/TN "^0.1. To interpret this result, there is not 
enough time for the a-process to accumulate a num- 
ber of seed elements and plenty of free neutrons are 
left even at the betyinninp; of the r-prncess. Conse- 



quflntly, the n/.s ra.tio becomes very high ~ IDfl. 



As for the neutron mass fraction, on the other 
hand, our value Xn = 0.159 is smaller than Woosley's 
model value Xn — 0.176 in trajectory 40 because 
low entropy is in favor of low neutron fraction. This 
may be a defect in our low entropy model. How- 
ever, the short dynamic time scale saves the situa- 
tion by regulating the excess of the seed elements as 
discussed above. These two effects compensate with 
each other to result in average mass number of seed 
nuclei {A) 95 and neutron-to-seed abundance ra- 
tio n/s ~ 100, which is ideal for the production of 
the third {A k. 195) abundance peak of the r-process 
elements in our model, as displayed in Fig. |l^. 

R-process elements have recent ly been detected in 
several metal-deficient halo stars ( ^neden ct al. 199^ ) 
and the relative abundance pattern for the elements 
between the second and the third peak proves to 
be very similar to that of the solar system r-process 
abundances. One of the possible and straightforward 
interpretations of this fact is that they were produced 
in narrow window of some limited physical condition 
in massive supernova explosions, as studied in the 
present paper. These massive stars have short lives 
~ 10^ yr and eject nucleosynthesis products into in- 



terstellar medium continuously from the early epoch 
of Galaxy evolution. It is not meaningless, there- 
fore, to discuss several features of our calculated result 
in comparison with the solar system r-process abun- 
dance distribution (Kappeler, Beer, & Wisshak 1989) 
in Fig. Although Kappeler et al. obtained these 
abundances as s-process subtractions from the ob- 



served meteoritic abundances (Anders and Grevesse 
T98|) for the mass region 63 < A < 209, the inferred 
yields and error bars for A = 206, 207, 208, and 209 
are subject to still uncertain s-process contribution. 
We did not show these heavy elements A — 206 — 209 
in Fig. |l^. 

Our single wind model reproduces observed abun- 
dance peaks around A k, 130 and A « 195 and the 
rare earth element region between these two peaks. 
However, there are several requirements to the wind 
model in order to better fit the details of the solar sys- 
tem r-process abundances in the mass region 120 A. 
The first unsatisfactory feature in our model calcula- 
tion is that the two peaks are shifted upward by 2 ^ 4 
mass unit, although overall positions and peaks are in 
good agreement with the solar system data. This is 
a common problem in all theoretical calculations of 



the r-proce ss nucleosynthesis ( Meyer et al. , Woosley 



et al. 19941) . The shift of the peak around A « 195 
is slightly larger than that around at A sa 130, which 
may be attributed to a strong neutron exposure as 
represented by n/s ~ 100 in our model calculation. 
The second feature is that the rare earth element re- 
gion shows broad abundance hill, but its peak position 
A Ki 165 in the data is not explained in our calcula- 
tion. It was pointed out by Surman et al. (1997) that 
the abundance structure in this mass region is sensi- 
tive to a subtle interplay of nuclear deformation and 
beta decay just prior to the freeze-out of the r-process. 
More careful studies of these nuclear effects and the 
dynamics of the r-process nucleosynthesis are desir- 
able. The third failure in the model calculation is the 
depletion around A « 120, which is also another se- 
rious problem encountered commonly by all previous 
theoretical calculations. This deficiency is thought to 
be made by too fast runaway of the neutron-capture 
reaction flow in this mass region. This is due to too 
strong shell effects of the = 82 neutron shell clo- 
sure, suggesting an incomplete nuclear mass extrap- 
olations to the nuclei with Z « 40 and A^ « 70 - 80 
which correspond to the depleted abundance mass re- 
gion A « 120. It is an interesting suggestion among 
many others (Woosley et al. 1994) that an artificial 
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smoothing of extrapolated zigzag structure of nuclear 
masses could fill the abundance dip around A « 120. 
This suggestion sheds light on the improvement of 
mass formula. 

Let us repeat it again that an overall success in 
the present r-process nucleosynthesis calculation, ex- 
cept for several unsatisfactory fine features mentioned 
above, is only for heavier mass elements 130 < A in- 
cluding the second {A w 130) and the third {A w 195) 
peaks. When one looks at disagreement of the abun- 
dance yields around the first {A « 80) peak, rel- 
ative to those at the third peak, between our cal- 
culated result and the solar system r-process abun- 
dances, it is clear that a single wind model is un- 
able to reproduce all three r-process abundance peaks. 
The first peak elements should be produced on differ- 
ent conditions with lower neutron-to-seed ratio and 
higher neutrino flux. I t has already bee n pointed 



out by several auth o rs ([Seeger et al. 1965|, [Kodama 



& Takahashi 1975| , [Hillebrandt, Takahashi, fc Ko- 
da 4ia 197^ that even the r-process nucleosynthe- 
sis needs different neutron exposures similarly to the 
s-process nucleosynthesis in order to understand the 
solar system r-process abundance distribution. In a 
single event of supernova explosion, there are several 
different hydrodynamic conditions in different mass 
shells of the neutrino-driven wind (Wooslcy ct al 



19c(|, [Witti, Janka, fc Takahashi 1994| ), which may 
produce the first peak elements. Different progenitor 
mass supernova or the event like an exploding accre- 
tion disk of neutron-star merger might contribute to 
the production of the r-process elements. Considera- 
tion of these possibilities is beyond our scope in the 
present paper. 

We did not include the effects of neutrino absorp- 
tion and scattering during the nucleosynthesis pro- 
cess in the present calculation. This is because these 
effects do not change drastically the final r-process 
yields as far as the dynamic expansion time scale Tdyn 
is very short. Using Eq. (^9|), we can estimate the 
specific collision time for neutrino-nucleus interaction 



0.082s -0.31s, 



(43) 



where the input parameters are set equal to L^^^i = 
10,61, = 15 MeV, and (a^) = 10-*^ cm^. Note that 
Tj/ « 0.082 s is the specific neutrino collision time 
at r = 52 km where the temperature of the wind be- 
comes T = 0.5 MeV at the beginning of the a-process, 
and T„ « 0.31 s for r = 101 km and T = 0.5/ e w 0.2 
MeV at the beginning of the r-process. These val- 



ues are larger than Tdyn — 0.0062 s that stands for 
the duration of the a-process by definition. There- 
fore, the neutrino process does not disturb the hydro- 
dynamic condition of the rapid expansion during the 
a-process. 

It is to be noted, however, that the neutrino pro- 
cess virtually makes the strong effect on the r-process 
for the winds of slow expansion. We have numer- 
ically examined Woosley's model (1994) of trajec- 
tory 40 to find Tdyn ~ 0.3 s. Meyer et al. (1998) 
also used Tdyn = 0.3 s in their simplified fiuid tra- 
jectory to see the neutrino-capture effects. This dy- 
namic time scale Tdyn ~ 0.3 s is larger than or com- 
parable to the specific neutrino collision time t^, in 
Eq. (p3[) . In such a slow expansion the neutrino ab- 
sorption by neutron (^) proceeds to make a new 
proton in the a-process. This proton is quickly in- 
terconvcrtcd into alpha-particle in the following re- 
action chain, p[n,^)d{n,^)i which is followed by 
t(p, n)'^IIe(ri, 7)^He and t(i, 2n)^IIe, and contributes 
to the production of seed elements. These radia- 
tive capture reactions and nuclear reactions are much 
faster than the weak process (|^) on newly produced 
proton from the process (H). The net effect of these 
neutrino processes, therefore, is to decrease the neu- 
tron number density and increase the seed abundance, 
which leads to extremely low n/s ratio. As a result, 
even the second abundance {A k, 130) peak of the r- 
process elements disappears, as reported in literature 
( [Meyer, McLaughlin, fc Fuller 1998| , [Meyer 1995| ). 
Details on the neutrino process will be reported else- 
where. 

We have assumed that electrons and positrons are 
fully relativistic throughout the nucleosynthesis pro- 
cess. However, the total entropy of the system may 
change at the temperature T <^ l/Snie where elec- 
trons and positrons tend to behave as non-relativistic 
particles. This might affect the nucleosynthesis al- 
though it does not affect significantly the dynamics 
near the protoneutron star. We should correct this 
assumption in the future papers. 

Finally, let us refer to massive neutron star. Large 
dispersion in heavy element abundance of halo stars 
has recently been observed. Ishimaru and Wanajo 
(1999) have shown in their galactic chemical evolu- 
tion model that if r-process nucleosynthesis occurs in 
either massive supernovae > 30Mq or small mass 
supernovae 8 — IOMq, where these masses are for 
the progenitors, the observed large dispersion can be 
well explained theoretically. In addition, SN1994W 
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and SN1997D are presumed to be due to 25 Mq — 
40 Mq massive progenitors because of very low '^^Ni 
abundance in the ejecta (^oUerman, Gumming, fc 



Lundqvist 1998, Turatto et al. 199S). These mas- 



sive supernova are known to have massive iron core 



> 1.8Mq and leave massive remnant (Turatto et al 
19f|8|). It is critical for the r-process nucleosynthesis 
whether the remnant is neutron star or black hole. 
Recent theoretical studies of the EOS of neutron star 
matter, which is based on relativistic mean field the- 
ory, set upper limit o f the neutron star mass at 2.2 
Mq (IShcn et al. 1998| ). 



5. Conclusion and Discussions 

We studied the general relativistic effects on neutrino- 
driven wind which is presumed to be the most promis- 
ing site for the r-process nucleosynthesis. We assumed 
the spherically symmetric and steady state flow of 
the wind. In solving the basic equations for relativis- 
tic fluid in Schwarzschild geometry, we did not take 
approximate method as adopted in several previous 
studies. We tried to extract generic properties of the 
wind in manners independent of supernova models or 
neutron-star cooling models. 

The general relativistic effects introduce several 
corrections to the equations of motion of the fluid and 
also to the formulae of neutrino heating rate due to 
the redshift and bending of neutrino trajectory. We 
found that these corrections increase entropy and de- 
crease dynamic time scale of the expanding neutrino- 
driven wind from those in the Newtonian case. The 
most important corrections among them proves to be 
the correction to the hydrodynamic equations. Both 
temperature and density of the relativistic wind de- 
crease more rapidly than the Newtonian wind as the 
distance increases without remarkable change of the 
velocity at r < 30km, where the neutrino-heating 
takes place efficiently. The lower the temperature and 
density are, the larger the net heating rate is. This 
is the main reason why the entropy in the relativistic 
case is larger than the Newtonian case. 

We also looked for suitable environmental condi- 
tion for the r-process nucleosynthesis in general rel- 
ativistic framework. We studied first the differences 
and similarities between relativistic and Newtonian 
winds in numerical calculations, and then tried to in- 
terpret their behavior by expressing gradients of the 
temperature, velocity and density of the system an- 
alytically under the reasonable approximations. We 



extensively studied the key quantities for the nucle- 
osynthesis, i.e. the entropy S and the dynamic time 
scale Tdyn of the expanding neutrino-driven wind, and 
their dependence on the protoneutron star mass, ra- 
dius, and neutrino luminosity. We found that more 
massive or equivalently more compact neutron star 
tends to produce explosive neutrino-driven wind of 
shorter dynamic time scale, which is completely differ- 
ent from the result of the previous studies in the New- 
tonian case which adopted approximation methods. 
We also found that the entropy becomes larger as the 
neutron star mass becomes larger. Since the larger 
luminosity makes the dynamic time scale shorter, the 
large neutrino luminosity is desirable as far as it is 
less than lO^^ergs/sec. If it exceeds 10^^ ergs/sec, 
only the mass outflow rate becomes large and the flow 
can not cool down to ~ 0.2MeV by the shock front 
r ~ 10,000km. As the result, the time scale becomes 
too long, which is not favorable to the r-process nu- 
cleosynthesis. 

Although we could not find a model which pro- 
duces very large entropy S ~ 400 as suggested by 
Woosley et al.(1994), it does not mean that the r- 
process does not occur in the neutrino-driven wind. 
We compared our results with Hoffman's condition 
and found that the short dynamic time scale Tdyn ~ 
6ms, with M = 2.0Mq and L^, = 10^^ ergs/sec, is 
one of the most preferable condition for producing r- 
process elements around the third peak {A ^ 195). In 
order to confirm this, we carried out numerical cal- 
culations of the r-process nucleosynthesis upon this 
condition by using fully implicit single network code 
which takes account of more than ~ 3000 isotopes and 
their associated nuclear reactions in large network. 
We found that the r-process elements around A ~ 195 
and even the heavier elements like thorium can be 
produced in this wind, although it has low entropy 
S 130. The short dynamic time scale Tdyn 6 
ms was found to play the role so that the few seed 
nuclei are produced with plenty of free neutrons left 
over at the beginning of the r-process. For this rea- 
son the resultant neutron-to-seed ratio, n/s ^ 100, 
in high enough even with low entropy and leads to 
appreciable production of r-process elements around 
the second(v4 « 130) and third {A « 195) abundance 
peaks and even the hill of rare earth elements between 
the peaks. 

Note that the energy release by the interconver- 
sion of nucleons into a-particles at T ^ 0.5 MeV pro- 
duces an additional entropy about AS ^ 14. This 
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was not included in our present calculation. We can 
make note that, taking account of this increase, the r- 
process could occur in the neutrino-driven wind from 
hot neutron star whose mass is smaller than 2.0Mq. 

One might think that short Tdyn brings deficiency 
of neutrino heating and that the wind may not blow. 
It is not true because the mass elements in the wind 
are heated by energetic neutrinos most efficiently at 
r <^ so/cm, while the expansion time scale Tdyn is the 
time for the temperature to decrease from T ^ 0.5 
MeV to 0.2 MeV at larger radii. The duration of 
time for the mass elements to reach 30km after leaving 
neutron star surface is longer than Tj^yn- There is 
enough time for the system to be heated by neutrinos 
even for Tiiyn as low as 6 ms. 

We did not include neutrino-capture reactions that 
may change Yg during the nucleosynthesis process. 
Since the initial electron fraction was taken to be 
relatively high Yg = 0.4 ^ 0.5, there is a possibil- 
ity that the final nucleosynthesis yields in neutrino- 
driven wind may be modified by the change in Y^. 
during the a- and r-processes. However, this is ex- 
pected to make a small modification in our present ex- 
pansion model with short dynamic time scale because 
the typical time scale of neutrino interaction is longer 
than Tdyn- We will report the details about the nucle- 
osynthesis calculation including neutrino-capture re- 
actions in forthcoming papers. 

It was found that the entropy decreases with in- 
creasing neutrino luminosity. This fact suggests that 
one cannot obtain large entropy by merely making 
the heating rate large. The cooling rate, on the 
other hand, does not depend on the neutrino lu- 
minosity. In the present studies we included two 
cooling mechanisms of the e~^e~ capture by free nu- 
cleons and the e~^e~ pair annihilation. As for the 
cooling rate due to the e'^e~ pair annihilation, only 
the contribution from pair-neutrino process is usu- 
ally taken into consideration, as in the present cal- 
culation. However, there are many other processes 
which can contribute to the total cooling rate. They 
are the photo-neutrino process, the plasma-neutrino 
process, the bremsstrahlung-neutrino process and the 



recombination-neutrino process (Itoh, Hayashi, & Nishikawa 
19f|5|). Indeed, if we double our adopted coohng rate 
artificially, we can obtain larger entropy. Details on 
the numerical studies of the cooling rate are reported 
elsewhere. The ra dial dependence of the heating rate 
is also important ( |Qian fc Wooslcy 1996| ) . Since both 
heating and cooling processes are critical to deter- 



mine the entropy, more investigation on the neutrino 
process is desirable. 

There are other effects which have not been in- 
cluded in the present study. They are, for example, 
the mass accretion onto the neutron star, the time 
variation of the neutrino luminosity, convection and 
mixing of materials, and rotation or other dynamic 
process which break spherical symmetry of the sys- 
tem. These probably important effects may make 
several modifications to the present result. However, 
we believe that our main conclusion that there is a 
possibility of finding the r-process nucleosynthesis in 
an environment of relatively small entropy and short 
dynamic time scale is still valid. We conclude that 
the neutrino-driven wind is a promising astrophysical 
site for the successful r-process nucleosynthesis. 
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Fig. 1. Outflow velocity (a) and temperature (b) in 
Schwarzschild geometry as a function of the distance r from 
the center of the neutron star for various mass outflow rate 
M, where neutron star mass M = 1.4Mq and neutrino lu- 
minosity Ly^, = 10 ''^ ergs/s are used. Long dashed curve is 
for the critical mass outflow rate Merit = 5.2681 X 10"*' Mq, 
in which the velocity becomes supersonic through the critical 
point. Fives curves denoted by 1 to 5 corresponds respectively 
to M = 5.25x10"'', 5.15x10-6, 5.0855x10-'', 5.0x10-", and 
4.8 X 10"® Mq. Calculated result denoted by "3" meets with 
our imposed boundary condition of T = 0.1 MeV at r = 10000 
km. Entropy per baryon S and dynamic timescale Tjyn, which 
correspond to each curves from 1 to 5, are tabulated in Table 1. 
Note that the temperature denoted by "5" does not decrease 
to T = 0.5/e MeV within 10000 km (see Table 1). 

Fig. 2. — The same as those in Fig.l, for the case of M = 
2.0Mq, Li^e = 10®^ ergs/s. Long dashed curve is for the critical 
mass outflow rate Merit = 1.2459 x IO-^Mq, and fives curves 
denoted by 1 to 5 correspond respectively to M = 1.245 X 
10""', 1.240 X 10-*, 1.225 x IQ-"*, 1.215 X IQ-"', and 1.195 X 
IQ- ^ Mq. Calculated result denoted by "1" meets with our 
imposed boundary condition of T = 0.1 MeV at r = 10000 
km. Entropy per baryon S and dynamic timescale r^yn, which 
correspond to each curves from 1 to 5, are tabulated in Table 1. 
Note that the temperature denoted by "5" docs not decrease 
to T = 0.5/c McV within 10000 km (sec Table 1). 

Fig. 3. — Outflow velocity u{r), in unit of lO^cm/sec, tem- 
perature T{r), in unit of O.lMcV, and baryon mass density 
Ph(x), in unit of lO^g/cni'^ , as function of the distance r from 
the center of the neutron star with the protoucutron star 
mass M = 1.4Mq, radius i? = 10 km, neutrino luminosity 
Lu = 10^^ ergs/s and initial density lO^^g/cm^. Solid and 
broken lines display the results in Schwarzschild and Newto- 
nian geometries, respectively. We choose the mass outflow 
rate M = 5.0855 x W~^Mq/s for the Schwarzschild case and 
M = 1.2690 X IO-^Mq for the Newtonian case. See text for 
details of the outer boundary condition on M. 

Fig. 4. — Entropy per baryon S(r) as a function of the dis- 
tance r from the canter of the neutron star. Solid and broken 
lines arc the same as those in Fig. 1 for the same set of the 
input parameters. 

Fig. 5. — Specific neutrino heating rate g(r) as a function 
of the distance r from the center of the neutron star, for the 
same set of the input parameters as those in Fig. 1. (a) To- 
tal net heating rate q. The solid and broken lines are for the 
Schwarzschild (denoted by general relativity) and Newtonian 
case, respectively, (b) Decomposition of the net heating rate 
into five different contributions from the heating process qi, 
53, and 55 (solid lines) and the cooling process 52 and 54 (bro- 
ken line) for the Schwarzschild case. See text for details of q;. 
(c) The same as those in (b) for the Newtonian case. 

Fig. 6. — Dynamic time scale T^yn (a) and entropy per baryon 
S (b) vs. neutron star mass M at 0.5MeV. Closed circles, con- 
nected by thick solid line, and open triangles, connected thin 
solid line, are the calculated results for the Schwarzschild and 
Newtonian case, respectively, by using the same set of the input 
parameters as in Fig. 1. Two broken lines are from Qian and 
Woosley (1996) in Newtonian case, which adopted an assump- 
tion of the radiation dominance (lower in T^yn, and upper in 



S) or the dominance of non-relativistic nucleon (upper in r^yn, 

and lower in S). 

Fig. 7. Dynamic time scale r^y^ (a) and entropy per baryon 
S (b) vs. neutrino luminosity Ly at T = 0.5 MeV. Thick and 
thin lines and two broken lines arc the same as those in Fig. 4. 
At the larger and of L,y ~ 10'''^ ergs/s, there is no solution 
to satisfy our imposed boundary condition, T = 0.1 MeV at 
r = 10000 km. See text details. 

Fig. 8. — Relation between entropy per baryon S and dy- 
namic time scale r^y^ for various combinations of the neutron 
star mass 1.2Mq < M < 2.0AfQ and the neutrino luminosity 
10®" < Lv < lO®'^ ergs/s. Solid and broken lines connect the 
same mass and luminosity. At the largest end of Li, ~ 10®^ 
ergs/s for each M, there is no solution to satisfy our imposed 
boundary condition, T = 0.1 MeV at r = 10000 km. Two zones 
indicated by shadows satisfy the approximate conditions, for 
Ye = 0.4, on which the successful r-process occurs (Hoffman et 
al. 1997) to make the second abundance peak around A = 130 
(lower) and the third abundance pealc around A = 195 (upper) . 
See text for details. 

Fig. 9. Seed abundances at T9 = 2.5 as a function of atomic 

immbcr A. Sec text for details. 

Fig. 10. Final r-proccss abundances (lines) as a function 
of atomic mass number A compared with the solar system r- 
process abundances (filled circles) from Kappeler, Beer, & Wis- 
shak (1989). The solar system r-process abundances are shown 
in arbitrary unit. Sec text for details. 

Fig. 11. The ratio of dynamic time scale r^yn to the time 
scale of typical a-proccss nuclear reaction t>j, rdyn/T^, vs. 
baryon mass density at T = 0.5 MeV, for various combina- 
tions of the neutron star mass 1.2Mq < M < 2.OM0 and the 
neutrino luminosity lO*'' < < 10®^ ergs/s. Solid and broken 
lines connect the same mass and luminosity. 
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Table 1 

Entropy and dynamic time scale for different M. Since the temperature in the 5th case for 

BOTH M = IAMq and M = 2.OM0 does NOT DECREASE TO T = 0.5/e MeV WITHIN 10000 KM, Tdyn IS NOT 

defined (see Fig. 1(b) and 2(b)). 







M 


entropy 


''dyn 






{IO-^Mq/scc) 


(fc) 


is) 


1.4M0,lO5iergs/s 


Merit 


5.2681 


116 


0.037171 




1 


5.2500 


117 


0.041304 




2 


5.1500 


120 


0.084335 




3 


5.0855 


123 


0.16455 




4 


5.0000 


126 


0.71569 




5 


4.8000 


135 




2.0MQ,10'^^ergs/s 


Merit 


1.2459 X 10^ 


138 


0.00507 




1 


1.2450 X 10^ 


138 


0.00618 




2 


1.2400 X 10^ 


139 


0.01088 




3 


1.2250 X 10^ 


141 


0.08962 




4 


1.2150 X 10^ 


143 


2.6272 




5 


1.1950 X 10^ 


146 
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